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We recently derived the Hamiltonian of fermionic composites by an exact procedure of bosoniza- 
tion. In the present paper expand this Hamiltonian in the inverse of the number of fermionic states 
in the composite wave function and give the necessary and sufficient conditions for the validity of 
such an expansion. We compare the results to the Random phase Approximation and the BCS 
theory and perform an illustrative application of the method. 
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The low energy physics of several many fermion sys- 
tems is determined by the formation of composites, ac- 
companied or not by their condensation. Examples of 
simple condensation (ground state built by a large num- 
ber of composite bosons in one single-particle state while 
all other occupation numbers are of order 1) can be found 
in superconducting materials, ultracold fermionic atoms 
in magnetic traps ultrasmall metallic grains [2] and 
atomic nuclei [31 H] • Examples of multiple condensation 
(ground state built by a large number of bosons in 2 
or more single-particle states) essentially occur in spa- 
tially separated systems as in the Josephson effect and 
superleak or in systems whose components differ by some 
internal quantum number of nongeometrical nature [S]. 
The theory of all such systems requires some bosonization 
procedure. For finite systems it is important to respect 
particle number conservation, and this becomes unavoid- 
able in the study of effects related to the total number 
of fcrmions, its parity and supersymmetry But re- 
specting this symmetry can be important also for infinite 
systems |S]. 

There is an infinite literature on this subject which we 
cannot review here, but to our kwnoledge there is no uni- 
fied mehtod which respects fermion number conservation 
and treats ground and excited states on the same foot- 
ing. We recently proposed an exact bosonization proce- 
dure which, as such, implements these features. Starting 
from the operator form of the partition function, after 
discretization of the Euclidean time we perform an in- 
dependent Bogoliubov transformation at each time slice. 
The time dependent parameters of the transformation are 
then associated with the fields of dynamical composite 
bosons in the presence of quasiparticles satisfying a com- 
positeness condition which avoids double-counting [7j. 
This formalism has been applied also to relativistic field 
theories, including gauge theories, at zero and finite tem- 
perature and fermion density [S]. 

In its first formulation [S] the method was restricted 
to complete bosonization (absence of unpaired fermions) . 
The resulting Hamiltonian was studied by an expansion 
in inverse powers of f2, the index of nilpotency of the 
composites, which is the number of fermionic states in 
their structure functions, assuming simple condensation 
with one-dimensional condensation sector. The conden- 



sation sector is the Fock subspace of composite bosons 
which contains the condensed boson and all other bosons 
whose mixing with the condensed one is allowed by sym- 
metries. We call free space the complementary subspace. 
In the present work we generalize the results of Ref . j9] 
to a condensation sector of arbitrary dimension. In this 
framework we can study BCS-BEC transitions in ultra- 
cold systems in magnetic traps [10] , effects related to the 
number of fermions and the parity of this number, shape 
transitions and shape coexistence jllj in atomic nuclei. 
One of the motivations of our work is indeed to derive the 
phenomenological Interacting Boson Model ^ of atomic 
nuclei starting from an effective nuclear Hamiltonian in 
a model space. Fermion Hamiltonians in a model space 
can be written in the form 
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cFkc. 
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c' , c are fermion creation-annihilation operators and 
we adopt a summation convention over repeated in- 
dices and a matrix notation, for instance cFkc = 



mim2 ■ 



e is the matrix of single- 
fermion energy, M the matrix of any interaction with 
external fields and Fk are the potential form factors 
normalized according to ti{Fl^_^FK2) = 2 il(5A'ii<'2 ■ The 
separable terms of the potential are assumed of particle- 
particle type: the particle-hole ones must be rearranged 
in this form. The Hamiltonian of composite bosons [9] 
is 

Hb =: tr{RB^eB) -^^gx^RB^ Fl^)tr{RFKB) 

K 

9Kti' \{R - 1)fI Fk - RB^fI^RFkB] : (2) 
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where the 

vn ^^^^ ' 



colons 



normal ordering, B 



mean 

R = (JL + B^B) and we omitted the inter- 
action with external fields Hj —: tr{RB^M B) : . In the 
above equations the Sk are destruction operators of com- 
posites with quantum numbers K and the matrices Bk 
are their structure functions. Unlike the composite op- 
erators Bk — c Bk c, they obey canonical commu- 
tation relations provided the conditions on the structure 
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functions reported below (4[5) are satisfied. In a sector 
of np fermions Hb must used with the constraint [9j 



tr (RB^'B) ^np. 



(3) 



We denote the quantum numbers of the composite bosons 
in the condensation/free sector by c, / fi respectively. We 
call "zero" the condensation mode and a the remaining 
states in the condensation sector. We will use the index 
K to represent both a- and /-bosons. 

The difficulty dealing with the bosonic Hamiltonian 
comes from the operator R which is a nonpolynomial 
function of bosonic creation-destruction operators. If the 
number of composite bosons is much smaller than we 
can overcome the difficulty by expanding R with respect 
to B^B. Otherwise a subtraction of the condensate is 
necessary before an expansion can be done [H]. Because 
only the "zero" mode will have a large occupation, we can 
make the number conserving approximation il^^ 6q6o ~ 
n/Q = i/Q, where n = \nF is the number of bosons. Now 
we can write the normalization conditions 



^trCoo = l,^tr 



= 5 



(4) 



where CKiKiKaKi... 



(5) 

.. and r 



^1 H- vBqBIj . They are derived as in Ref.j9j and coin- 
cide with those reported there if the condensation sector 
is one-dimensional. The subtraction parameter h' (a pri- 
ori different from vq) is determined by the constraint ([3| 
which for simple condensation using Eqs. (4j5l becomes 
n{vQ — v) = terms of order 17°. This equation is satisfied 
hy V = vq + terms of order Q,"^ which give contributions 
of order 1 which depend on the number of fermions but 
not on their state. Following Ref.[5] Hq can be expanded 
in powers of fi^i. Neglecting terms of order f2~2 it be- 
comes 



Hb = 

H 

where 



£oo n + £o^ blb„ + £„o blbo + S^Ji^ ^TiJ^2 



B.C. 
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We used the definitions 

TK1K2...K21-1K2 



K2 



.c 
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tr [CkiK2---Ck21-3K21- 



xTBk2,_,F},^^ 
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where I is any integer. The above equations are valid 
for an arbitrary choice of e as a consequence of the con- 
straint ([3]) on the number of fermions. All the constants 
appearing in Hb are of order Then assuming n ~ 
the first term is of order fi, the second and third are 
of order fis (because the matrix elements of the opera- 
tors bg b„ are of order ^Jn) , and the other terms are of 
order fi^. Estimating the order of the different terms 
we assume that the coupling constants scale according 
to gii ^ ri~^(5e, where Se is the single fermion energy 
spreading. This assumption is necessary to get finite en- 
ergies in infinite systems. Two comments are in order: 
i) all the constants appearing in Hb vanish unless both 
quantum numbers K1K2 belong either to the condensed 
or to the free sector because of different symmetries in 
these sectors, ii) the zero-boson has a strong mixing with 
the cr-bosons due to the terms of order fis . But if it is 
the true condensed boson such mixing must be absent 



£o<j — £a-o — . 



(8) 



As we will see this decoupling prevents the occurrence 
of boson-particle-boson-hole states in the ground state, 
and it is therefore similar to an ordinary selfconsistency 
condition. 

We can now derive a number of results without numer- 
ical calculations by introducing phonon operators. This 
can be done in two steps. First we introduce the phonon 
operators A-^ ~ b-^ &q and rewrite the Hamiltonian 
accordingly 



Hb 



£0 + Soon + Sj^^j^^aL Aj^ 
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K1K2 Ki K2 



H.C 



(9) 



This is the Bogoliubov Hamiltonian for a superfluid 
bosonic system and can be diagonalized by a bosonic Bo- 
goHubov transformation = <|t(7t _|_ |>]/t ^ ^ = (j^ 

This transformation becomes simple when, due 
to symmetries, there is a one-to-one correspondence be- 
tween the quantum numbers K2 and Ki: K2 — K2{Ki). 
This is guaranteed only for the free sector, but we as- 
sume it to be true for all noncondensed modes in or- 
der to illustrate some general features. Then Sj^^j^^ — 



K1K2 

mation Hg takes its final form 



^Ki^K2{Ki),Ki- 



After transfor- 



HB=£o + J2y^-^K+^oon + Ej^^\^^j^ (10) 



K 



where Ej^ = y^f ^ — Vq VV^. The ground state is the 

vacuum of the phononic operators <i>, which written in 
terms of the bosonic composites is 



\GS) 



exp 



1 



2n 



{ViU* 



I GS) is a coherent state of phonon pairs which are 2p-2h 
boson states and therefore 4p-4h fermion states built on 
the 60 condensate, at variance TT with the RPA. The 
above results have been derived under the condition that 
the occupation of noncondensed bosons be of order 1 
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K 



«n. 



(11) 

This is the necessary and sufficient condition for the va- 
lidity of the assumption of simple condensation and of 
the expansion. So this is also the condition for con- 
densation in a finite system. 

The structure functions Bx must be determined by 
a variational calculation on ground state and phonon 
energies under the normalization and decoupling condi- 
tions. This can be done in the following way. We write 
Bq = Bq + ^iSBq. We solve the variational equation 



for Bq neglecting the first 2 terms of (10) which are of 



order 1 and use the result to determine the correction 
SBq. In the present approximation in which we neglect 
terms of order f2~ 2 this correction will alter only the con- 
tribution of the term Soon. A comparison with the BCS 
theory can clarify the meaning of this procedure. The pa- 
rameter £00 is indeed the BCS energy per particle in the 
form of the quasichemical equilibrium theory [13 . To see 
this we must sketch an essential point in the derivation 
of Hb. The time-dependent Bogoliubov transformation 
changes the fermion vacuum into the quasiparticle vac- 
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uum \B) — exp (^bxBl^^ |0), where the Bk = cB 

are composite operators and the bx time-dependent holo- 
morphic variables. The partition function involves an in- 
tegral over these variables which acquire the meaning of 



dynamical bosonic fields with fermion number 2, so that 
all the terms in the expansion of \B) have fermion number 
zero. A further transformation of the partition function 
into a trace over a bosonic Fock space leads to the ex- 
pression ([2|. The coherent states \B) can be identified 
with the BCS states setting 



(Bo)r 



n 



(Bo) 



2 



with standard BCS notation. But if one is ready to 
break fermion number conservation one can set 6g « n. 
Then the BCS normalization + v"^ = 1 is identi- 
cally enforced and the first of the normalization condi- 
tions ffl becomes the BCS condition on fermion number 

One can easily check that the variation of £00 with re- 
spect to -Bo yields the standard gap equation. But this 
equation must now be solved under the decoupling con- 
ditions ([8|. These conditions are obviously absent when 
the condensation sector is onedimcnsional, in which case 
one can use results obtained in the BCS theory to eval- 
uate the phonon energies and 5Bq which corresponds to 
the contribution arising from projection of the BCS wave 
function (but it includes the contribution of states in the 
free sector). We now show how our method can be ap- 
plied in practice to the study of shape transitions in nu- 
clei and the BEC-BCS transition at zero temperature, 
which can be regarded as the transition from molecu- 
lar to Cooper structure of the condensed boson. In fact 
the convenience of the quasichemical equilibrium theory 
in this context has been advocated long ago [2]. In 
the principal frame of inertia of the nucleus the zero- 
mode will be a superposition of &iM-bosons with M = 0, 
the (7-modes will be all the other independent superpo- 
sitions of modes with M — Q and the free modes will 
have M ^ 0. To determine the structure functions we 
must minimize the ground state and the phonon energies 
under the decoupling and normalization conditions ( 8|4 
[5]). The problem simplifies considerably if we assume 
the condensation sector to be 2-dimensional, including 
only the states L = 0,2, the s- and d- bosons of the 
Interacting boson model. As a further simplification 
we assume the nucleons to live in a single j-shell and 
have constant single-particle energy and monopole and 
quadrupole paring interactions. The energy is therefore a 
function of the ratio of couplings ^ = 52 /5o- In this model 
space 17 = J + I and the potential form factors are pro- 
portional to Clcbsh-Gordan coefficients {FLM)m-im2 ~^ 
jmijm2\LM > yjlj + 1 . Then the structure matrices in 
the condensation sector can be parametrized according 
to Bq = 70F00 -t- 72^20 , ^<T = A^(cos/3i^oo + sin/3F2o), 
and all the parameters are determined by the constraints. 
/3 and A/", are determined by Eqs. [i] Eqs. ^ and ([s]) 
give rise respectively to a closed and open contour lines in 
the 70,72 plane as shown in Fig.l where only the upper 
part is shown (the lower one can be obtained by a parity 
inversion). The intersection for which £00 is minimum 
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FIG. 1: Contour plots of the first of normalization conditions 
(|4|(full circles) and of the decoupling condition (|8| (empty 
circles) for ^ = 1.5 and vo = 0.5. The minimum energy is for 
a spherical shape, at 72 = 0,70 = 1.4. 
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FIG. 2: The critical line separating the deformed phase (up- 
per part) from the spherical one (lower part). 



determines the values of the 7-parameters and therefore 
the shape of the nucleus. In Fig. 2 we plot the criti- 
cal line which separates the deformed from the spherical 



shape. When i/q approaches 1 the shell is fully occupied, 
and because of the Pauli principle the bosonic approx- 
imation cannot be valid. However we have a tentative 
interpretation of the steep rise of the critical line in this 
region. For vq = ^ there is a unique nuclear state which 
is spherical. Therefore approaching 1, larger and larger 
values of the strength of the quadrupole pairing are nec- 
essary to get a deformed shape. Below the critical line the 
ground state is a condensate of bosons with zero angular 
momentum (seniority limit). The condensation sector is 
onc-dimcnsional and 



where 



Er 



1 



(12) 



(13) 



Note that the Hamiltonian is expressed in terms of the A- 
phonons, because the coefficients W vanish. The phonon 
energies are even functions of the deviation from half 
filling and vanish for ^ = 1, close to the critical line. 
For ^ = we reproduce the spectrum of the pairing 
model [12]. Inclusion of a nonconstant single-particle 
energy does not require any modification of the proce- 
dure outlined. Instead if we increase the model space 
and/or the condensation sector, the structure functions 
will depend on more parameters and we will have to solve 
a constrained gap equation. The resulting Hamiltonian 
can be directly compared to that of the Interacting Boson 
model in the intrinsic frame. 
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